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I present results of a numerical calculation of the effects of light quark- 
^\ . antiquark pairs on the linear heavy-quark potential in light-cone quantized 

On ' two-dimensional QCD. I extract the potential from the QQ component of 

the ground-state wavefunction, and observe string breaking at the heavy- 
^s ■ light meson pair threshold. I briefly comment on the states responsible 

I ' for the breaking. 
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1 Introduction 

Pure Quantum Chromodynamics exhibits linear confinenient between charges at 
large distances, binding them with a gluonic string, as observed some time ago 
in lattice simulations of the potential between static quarks |l|] . Stretching the 
string sufficiently in the presence of light quarks should produce quark-antiquark 
pairs, breaking the string and screening the charges; attempts to clearly demon- 
strate this effect in lattice QCD are ongoing ||]. The purpose for this study is to 
determine whether string breaking can be observed in a relatively simple model, 
two-dimensional QCD [0, using a very different numerical scheme: solving for 
eigenstates of the light-cone (LC) Hamiltonian for the system confined to a box 

i- 

Early interest in LC quantization focused on high-energy scattering, due in 
part to the boost invariance of LC wavefunctions and their direct connection to 
parton distribution functions jsl. However, the suggestion of a simplified vac- 
uum and successes with two-dimensional models generated some recent interest 
in using this scheme to perform nonperturbative calculations of static properties 
in QCD [pj. In addition to their connection to parton distributions, LC wave- 
functions have a well-defined heavy-quark limit wj and a simple connection to 
quantum mechanical wavefunctions, which I will make use of in extracting the 
heavy-quark potential below. 

In this formalism, one initializes fields and their commutators at equal LC 
time, x~^ = x'^ + x^ = 0. The Hamiltonian P~ = P^ — P^ evolves states in a;+, 
while the LC momentum _P+ = P'^ + P'^ is kinematic and conserved in interac- 
tions. Diagonalizing P^ in the space of states defined by acting with creation 
operators on the vacuum produces the masses and wavefunctions for hadrons 
and multiple-hadron states. All hadronic matrix elements may be expressed as 
integrals over these wavefunctions. 

In order to accomplish this numerically, I confine the system to a box of 
length L in x'^ and impose antiperiodic boundary conditions on the fields. As 
a result, each quark carries fc+ in odd half integer units 27r/L. Because all 
fc"*" are positive and must sum to the total P"*", the discrete momentum K = 
{L/2tt)P^ determines both the resolution of the wavefunctions and also controls 
the number of quanta possible in a state. K —^ oo defines the continuum limit. 

The program which simulates this system allows for any number of colors 
and flavors. Once given conserved quantities such as baryon number, total 
momentum, and individual flavor numbers, it generates all states consistent with 
these, though it allows for additional restrictions. It computes the Hamiltonian 
in this basis and diagonalizes it, generating the complete set of masses and 
wavefunctions. That it produces all states is both an advantage and curse, 
as the number of states grows roughly exponentially with K. See J. Hiller's 
contribution to these proceedings [p| and references therein for a discussion of 
an intelligent way to focus on the lowest lying states. 

While fully relativistic, the use of a Hamiltonian, a Fock space description 
of states, and wavefunctions makes it similar in form to nonrelativistic many- 
particle quantum mechanics. A final similarity is the ability to work with a 



simple vacuum. Roughly, because all particles carry positive fc"*", it is not pos- 
sible for their momentum to sum to the zero P"*" of the vacuum; therefore, the 
vacuum has no particles. The true situation is more involved, and whether this 
argument holds depends on what is being computed. Extracting vacuum prop- 
erties in particular can be subtle B. But when applicable it is an enormous 
simplification. 

One of the advantages of LC quantization absent in the older infinite mo- 
mentum frame formulation is that it selects a quantization surface rather than 
a specific reference frame. As a result, it has a well-defined nonrelativistic limit 
for heavy quarks, in which LC wavefunctions go smoothly into the usual wave- 
functions of nonrelativistic quantum mechanics. Because I wish to extract the 
potential felt by two nearly static quarks, I include only one heavy flavor. Heavy 
in two-dimensional QCD means that TOq is large compared to the dimensionful 
coupling g; in four dimensions, the comparison would be to Aqcd- 

To see the Schrodinger equation these satisfy, consider the projection of the 
LC equation 

M^ 10) = P+P- 10} (1) 

onto the \QQ) subspace. States in this subspace have the form 

dx cj,,,^{x) h^\xP+) dt((l - x)P+) |0) (2) 



where x gives the fraction of P+ carried by the quark with probability \4'qq{x)\'^ i 
and the color index c contracts to form a singlet. While the complete wavefunc- 
tion contains states with extra QQ pairs, these are suppressed by factors of 
1/?tiq. The result is 't Hooft's equation [g|: 
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Here g^ = g'^{N'^ — 1)/2N, and the principal value prescription defines the 
integral. Though restricted to the two quark sector, this equation is still fully 
relativistic; in fact, it describes mesons to leading order in large N. 
In the large-rriq limit, the momentum fraction carried by the quark. 
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to lowest order in the relative momentum q, and the LC wavefunction 0qq(x) 
becomes sharply peaked around 1/2. Introducing a nonrelativistic momentum- 
space wavefunction 

V'(g) = 0QQ(l/2 + <z/2mQ) (5) 

focuses on deviations from this peak. Defining the nonrelativistic energy E = 
M — 2?71q and reduced mass fi = ?71q/2, extending q from irriQ to infinity. 



and Fourier transforming tp{q) to the relative position r turns Eq. (pf) into the 
Schrodinger equation mu, UU H 
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V'(r) = Ei;{r) . (6) 
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The relation 

^[e-'^'--l] = -^|r|. (7) 

converted the integral over q into a linear potential in position with string tension 
g^/2. The solutions are Airy functions; see Ref. jl^l for details. 

The computer code produces masses and LC momentum-space wavefunc- 
tions on whose Fourier transforms Eq. ffl) holds in the large- ttjq limit. There- 
fore, computing %lj(r) and d'^ip(r) by transforming ip{q) and q'^ip{q) respectively 
and defining Q 

_ 2fiE ^P{r) + dl jjjr) 

should reproduce a linear potential in that limit. In particular, the degree to 
which Vqq conforms to |r| indicates how accurately a nonrelativistic Schrodinger 
equation with a static potential describes this meson. Figures |l| and y display 
potentials extracted in this manner for two values oi g^/niQ. These reconstruct 
g'^\r\/2 fairly reliably at the relatively modest momentum K = 15, although the 
effects of discretization appear in the turnover of Vqq near its edges. 

2 String Breaking 

In one spatial dimension, there are no transverse directions in which the electric 
field may spread. As a result, the Coulomb potential is linear; in A~^ = gauge, 
it is all that remains of the gauge field. Because the string is present in the 
classical Hamiltonian, all the work on the quantum level goes into breaking it, 
making it a convenient laboratory to study this effect. To do so, I introduce 
an additional light-quark flavor and study the effect the presence of extra light- 
quark pairs has on the potential defined in Eq. (0). More specifically, because 
the program includes all states consistent with quantum numbers specified, I 
include three flavors: two different heavy quarks with identical mass and one 
light quark, but I fix their respective flavor numbers to 1, —1 and 0. This ensures 
that every state contains at least one heavy QQ pair. 

Figures |l| and y display preliminary results for the potential in the presence 
of light quarks at couplings of (^n/^^q — -164 and .247 and a light quark mass 
of TOq = .001/toq. The total discrete momentum K = 15 for these, leading to 
15 points in the QQ wavefunction. I used an invariant mass cutoff to exclude 
irrelevant extra heavy-quark pairs, and also restricted by hand the state space 
to a single additional qq pair in order to keep the number of states manageable; 
typically, about 1300. Testing this truncation by allowing additional pairs at 
lower K had little effect, suggesting that the production of a single light-quark 



pair was predominantly responsible for breaking the string. The data displayed 
came from a run of roughly half an hour on a DEC alpha 3000/700 workstation. 
Along with the ground states from which these potentials were extracted, the 
program produced the entire set of roughly 1300 masses and wavefunctions for 
this system. The LC Hamiltonian has a relatively simple structure, with the 
coupling and masses appearing as overall factors before the free and interacting 
terms. Almost all of the computational effort goes into producing a complete 
orthonormal set of states and evaluating the Hamiltonian matrix. Therefore, by 
storing separately the free and interaction matrices, reproducing the spectrum 
for different couplings and masses is essentially free. 

Both the linear potential for the pure QQ system and its breaking in the 
presence of light quarks are clearly evident. The horizontal line in each plot 
indicates the additional cost in energy to produce a pair of heavy-light mesons, 
showing that the string breaks about where expected. To obtain these, I com- 
puted the mass of the lightest Qq mesons with the same parameters but half the 
value of K. This approximates that carried by each heavy-light meson produced 
in the QQqq system. 

The ground state includes in its wavefunction not only the \QQ) states from 
which V"qq is extracted, but also higher- Fock states \QQqq) with additional 
light quark pairs. By examining their magnitude, it is possible to gain some 
insight into which states are most responsible for breaking the string. For the 
gN/rriQ = .247 case, states in which each light quark forms a singlet with a heavy 
quark and the LC momentum is carried almost entirely and equally by the heavy 
quarks are predominant. This is in accord with the tendency for Fock-state 
constituents to move with equal velocity ||l^ . These account for approximately 
.5% of the total probability. Figure]^ shows the difference between the position- 
space wavefunctions for g^/mQ — .247, before and after the inclusion of light 
quarks. These exhibit the spread away from r = which results from charge 
screening by the light quarks. 

In these figures, I computed Vqq from the ground state wavefunctions. Note, 
however, that Eq. (H) applies equally to excited states. Vqq extracted from these 
is consistent but significantly degraded. 

There are some tradeoffs evident in these preliminary results. The larger 
the coupling, the more dramatic the breaking and longer the extent of the 
wavefunction. But the larger coupling leads to a less relativistic QQ system, 
and a wider wavefunction in momentum, making it more susceptible to finite 
box size errors. It also increases the importance of additional qq pairs thus far 
neglected. 

3 Future Work 

In addition to simply increasing K, there are a number of ways of improving 
the accuracy of this study. Fitting the wavefunctions to physically motivated 
continuous functions would ameliorate some of the discretization and finite-size 
errors. The range of Vqq, which is limited by this method to the width of the 
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Figure 1: The potential defined in Eq. (@) at g^/rriQ — .164 for two heavy 
quarks (O), and after the inclusion of one light flavor with m^/niQ — .001 (+). 
The discrete momentum K — 15 for both cases. The horizontal line gives the 
additional energy 2MQ^ — 2mQ required to produce two heavy- light mesons, each 
computed at i^ = 8. The square root of the string tension gl/2 fixes the units 
for Vqq and r. 



wavefunction, could also be extended with the use of static sources, though 
this would mean abandoning the momentum conservation currently built into 
the code. Finally, an intelligent method for importance sampling of states, and 
perhaps the use of effective interactions to replace states excluded based on 
energy, would certainly improve the efficiency of this program. 
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Figure 2: The potential defined in Eq. (g) at gi^/mQ — .247 for two heavy 
quarks (O), and after the inclusion of one light flavor with m^/niQ ~ .001 (-I-). 
The discrete momentum K — 15 for both cases. The horizontal line gives the 
additional energy 2MQ^ — 2mQ required to produce two heavy- light mesons, each 
computed a,t K — 8. The square root of the string tension g^/2 fixes the units 
for Vqq and r. 
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